Measuring and engineering entropy and spin squeezing in weakly linked Bose-Einstein 

condensates 
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We propose a method to infer the single-particle entropy of bosonic atoms in an optical lattice 
and to study the local evolution of entropy, spin squeezing, and entropic inequalities for entangle- 
ment detection in such systems. This method is based on experimentally feasible measurements 
of non-nearest-neighbour coherences. We study a specific example of dynamically controlling atom 
tunneling between selected sites and show that this could potentially also improve the metrologically 
relevant spin squeezing. 
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The quest for novel cooling schemes to control the en- 
tropy of ultracold atoms in optical lattices is attracting 
considerable interest because thermal and quantum fluc- 
tuations limit the use of these systems for quantum sim- 
ulation or quantum metrology. For example, the experi- 
mental observation of magnetic ordering, a milestone for 
the quantum simulation of spin systems, is hindered by 
finite entropy in the system In metrology, finite 

entropy limits the amount of achievable spin squeezing, 
a useful resource for quantum-enhanced high precision 
measurements Not only controlling the entropy is 

challenging; it is also difficult to measure on a micro- 
scopic level. Only in the low atomic filling regime of the 
Mott-insulating phase a mapping of observable on-site 
atom number fluctuations to single-particle entropy has 
been achieved [10|, stimulating interest in entanglement 
detection by entropy measurements Recent ex- 

periments show that individual atoms can now be manip- 
ulated on a single-spin level at individual lattice sites [13l| 
and, in principle, the entropy can be locally engineered 
in order to effectively cool the system 14 1. 



We propose techniques for accurate local detection and 
control of a single-particle entropy of a bosonic atomic 
gas in an optical lattice in the large filling regime, opening 
up novel avenues for estimating entropy-based quantities 
and engineering entropy. In the case of non-negligible 
tunneling of atoms between adjacent sites, access to the 
intcr-site coherence is crucial for the determination of the 
entropy. An essential ingredient of our proposed entropy 
measurement is a technique for a spatially resolved mea- 
surement of long-range coherence of the lattice system 
based on a matter-wave homodyne measurement with 
respect to a reference condensate. Using this method 
we theoretically show how the entropy and entropic in- 
equalities of entanglement detection can be accurately 
estimated in experimentally realistic cases 0, m by the 
atom number in individual sites and the relative phase 
coherence between the atoms in different sites. The pro- 
posed coherence measurement may also lead to improved 



detection of spin squeezing of atoms between different 
sites, with potential applications to high-precision mea- 
surements. It also allows the detection of spin squeezing 
between non-nearest-neighbour sites and other spatially 
separated regions. Moreover, we study a specific exam- 
ple of controlling entropy by locally tailoring the trapping 
potential. We show how the coupling between the high- 
entropy regions and the rest of the system is adjusted by 
suddenly applying laser barriers between the central and 
outermost sites, altering the entropy distribution and af- 
fecting spin squeezing in the system. An illustration of 
the proposed scheme is shown in figure [T] together with 
the variation of the entropy and the on-site atom number 
fluctuations in the initial state: The barriers suppress the 
interactions between the atoms in the central sites and 
the outer-well atoms that exhibit stronger thermal and 
quantum fluctuations. We present example simulations 
in which injecting the laser barriers leads to reduced spin 
fluctuations and improved spin squeezing. We also find 
spin squeezing between atoms occupying spatially distant 
sites. 

We consider Bose-condensed atoms confined in an elon- 
gated trap in which we neglect density fluctuations in the 
radial direction. In the axial direction the atoms expe- 
rience a combined harmonic and optical lattice potential 
of a few sites, forming an array of weakly linked Bose- 
Einstein condensates (BECs). As the lattice height is 
increased the atom number fluctuations in each site are 
reduced and the system can exhibit metrologically rele- 
vant spin squeezing between the atoms in adjacent sites 
. Each lattice site occupies a multi-mode BEG and the 
atom number statistics in such a system is influenced by 
interactions between several modes in each site [T^ (in 
the simulations we include nine vibrational modes per 
site). By numerically solving for the amplitudes of the 
wavefunctions \faj), determined by the lattice sites a 
and vibrational levels j, we can construct the atom num- 
ber fluctuations, phase coherence, and the entire density 
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FIG. 1. (a) A schematic illustration of the relative phase 
measurement scheme, (b) Entropy engineering in the optical 
lattice where the tunneling of the atoms from the outermost 
sites to the rest of the system is suppressed, e.g., by apply- 
ing symmetrically-positioned, tightly-focused, blue-detuned 
laser beams. We show the combined lattice potential and 
the optical barriers (black solid and dashed lines), the initial 
atom density before the sites are decoupled (solid grey line), 
the single-particle entropy Sa/{na) per atom in site a (green 
dashed line with stars) , and the normalized on-site atom num- 
ber fluctuations (Ano,)^ /(ria) (red solid line with circles) in 
the initial state for each lattice site. The temperature is as- 
sumed to be 4 uK. 

matrix of the system, 

/5= ^Paj,^i|¥'aj)(¥',9,i|- (1) 

This can then be used to evaluate the single-particle von 
Neumann entropy 

5 = -Tr(plogp). (2) 

The density matrix elements Paj,pi are determined by 
the mode populations and their relative phase coherence 
(see p) . 

Experimentally, it is challenging to measure popula- 
tions of the individual vibrational levels or the relative 
phase coherence between them. In order to circum- 
vent the need to gather such detailed information, we 
will show that in experimentally realistic situations the 
single-particle entropy may be estimated by the atom 
number and the relative phase coherence between the 
atoms in different sites that are obtained by averaging 
over the vibrational level structure in each site. 



In the experiments [l^ a high-precision optical ab- 
sorption imaging provided site-resolved detection by in- 
tegration of the imaged atom density. Moreover, local 
interference measurements were performed after a short 
condensate expansion time allowing only the atoms from 
adjacent sites to overlap. The relative phase coherence of 
the atoms between the adjacent sites was then inferred 
from the phase variance of the interference pattern. In 
order to detect the long-range phase coherence in the 
lattice and extract sufficiently accurate information of 
the entropy, we propose a matter-wave homodyne mea- 
surement scheme for the atoms: We consider a system 
where the bosonic atoms in the lattice are surrounded 
by a BEC in a different internal state, e.g., a different 
hyperfine level for alkali-metal atoms or a metastable 
triplet state for alkaline-earth or rare-earth metal atoms 
(figure [T]). The atoms in this second internal state are 
assumed to experience a much weaker lattice potential 
or weaker interactions, so that their phase coherence is 
well preserved over the entire lattice length. Experimen- 
tally, inter-species interaction might be controlled by a 
Feshbach resonance or by adjusting the spatial overlap 
between the two species, e.g., with superlattices (l6l.[l7[. 
The surrounding BEC serves as a common phase refer- 
ence in analogy to local oscillators in quantum optical 
homodyne measurements. The phase coherence of the 
atoms in each lattice site can be locally determined by 
interfering the atoms with the reference condensate, e.g., 
by Raman transitions between the hyperfine levels [18|. 
Due to negligible phase fluctuations in the reference con- 
densate, local interference measurements of the lattice 
atoms with the reference condensate provide information 
about the relative phase fluctuations between the atoms 
in distant sites. The proposed scheme has the advan- 
tage of precise coherence readout by particle counting as 
experimentally problematic shot-to-shot fluctuations of 
the magnetic field (which lead to excess phase fluctua- 
tions) do not disturb such measurements as long as they 
are spatially homogeneous. The measurement method 
could also be suitable for 2D lattice systems when com- 
bined with recently developed high resolution imaging 
techniques [lol [l9|. 

In order to demonstrate how the long-range coher- 
ence and atom number detection can be used to infer 
locally a single-particle entropy and spin squeezing in an 
experimentally realistic system we study a specific ex- 
ample of dynamically adjusting the atom tunneling be- 
tween central and outermost sites (see figure [T|). We 
assume that the atoms are initially confined in a ther- 
mal equilibrium state in a shallow lattice. We then sim- 
ulate the resulting dynamics when the coupling of the 
outermost sites to the rest of the system is suppressed 
by a rapid injection of laser barriers, followed by a slow 
ramp up of the lattice potential. The increase in the lat- 
tice depth results in reduced atom number fluctuations 
and stronger spin squeezing between atoms in adjacent 
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sites. The laser barriers alter the entropy distribution 
in the system, as the initial thermal fluctuations in the 
trap are not uniform (figure [1]). Our numerical simula- 
tions are based on the truncated Wigner approximation 

(TWA) [IB-lii], using an approach similar to the one in- 
troduced in Ref. [15[. Here thermal and quantum fluc- 
tuations of the atoms in the stochastic initial state are 
calculated by self-consistently solving the ground-state 
and excited-state populations within the Hartree-Fock- 
Bogoliubov approximation [i^, S| • During the time dy- 
namics the field amplitudes in the wavefunction basis 
are obtained by projecting from the numerically 
calculated stochastic Wigner field (seCp). 

We take the experimental parameters that were used 
to observe spin squeezing between the atoms in adja- 
cent sites in which case N = 5300 atoms were con- 
fined in a combined ID lattice potential, with the spacing 
d ~ 5.7/zm, and an elongated harmonic trap with the fre- 
quencies a; ~ 27rx 21 Hz and wj^ ~ 27rx427 Hz {lu <^ ujj_), 
so that about 95% of the atoms occupied the six central 
sites. In Ref. [l^ the ID TWA model provided a good 
qualitative agreement with the experimental findings of 
the on-site and the relative atom number fluctuations. 
Here we use the same approach with the potential 
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The initial lattice height of 24 (with the recoil energy 
Er — h^ir'^ /2m(P) is slowly turned up to 72 Er. We 
study two different ramping speeds 15.6 Hz/ms and 17.2 
Hz/ms. The strength of the nonlinear atom-atom inter- 
action is given by gmN = iSlhiul, where gio — 2huj±a, 
a is the s-wave scattering length, and / = {h/mcoY^^. 
We consider the selection of the four central sites by in- 
jecting narrow, blue-detuned laser beams, or a scissors 
potential, before the ramping up of the lattice. The laser 
potential can be modeled by two symmetric Gaussian in- 
tensity distributions centered at ±Xb 



Vb = SbEn, 



exp 



{X - Xbf 



exp 



(x -I- XbY 



(4) 

with the waist 2db — 600nm (1/e^ intensity radius) that 
can be achieved by diffraction-limited focusing of a laser 
with wavelength in the visible light region. We vary the 
barrier height Sb and the cutting is fast compared to any 
other time scale so that the system has no time to relax 
during the process. 

We evaluate the entropy per atom Si in the central four 
wells and and compare this to the entropy per atom St 
in the full system. The effect of different strengths of the 
scissors potential is shown in flgure[2]by displaying the 
entropy in different cases at the end of the lattice ramp- 
ing. The simulation results of the barrier injection pro- 
vide an example that demonstrates how a single-particle 
entropy can be engineered in a nonequilibrium process by 
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FIG. 2. Estimates (dashed lines) derived from ((5| and exact 
results (solid lines) for the total entropy per atom St (lines 
with crosses) and internal entropy per atom Si (lines with 
circles) for different strengths of the scissors. Top row: Initial 
temperature T — 4nK for the ramping speeds (from left to 
right) 15.6 Hz/ms and 17.2 Hz/ms. Bottom row: The same 
for T — 0. The curve for the total entropy per atom is shown 
together with the sampling error calculated numerically over 
800 individual stochastic realizations (see,^. The insets show 
the ratio Si /St (exact results) for each case. 



laser barriers. Without scissors Si <^ S^^ since in a com- 
bined optical lattice and harmonic trap thermal phonon 
excitations and entropy are mostly concentrated on the 
outermost sites. Decoupling those sites from the rest of 
the system therefore reduces the effect of higher excita- 
tion modes. For weak scissors (s{, < 60), the entropy does 
not notably increase and 5*1 remains lower than the value 
of St without the potential barriers. In the insets in fig- 
ure [H we display the ratio Si /St- We find an increasing 
total entropy (but decreasing Si/St) for strong scissor 
potentials owing to excitations induced by injected bar- 
riers. A detailed investigation of these excitations and 
resulting entropy waves, which is beyond the scope of 
the present work, could in itself provide an interesting 
further study of entropy phenomena in a coupled multi- 
mode EEC system. 

As it is not practical to measure the entropy by de- 
tecting all the vibrational mode amplitudes, we propose 
an entropy estimate based on the lattice site occupation 
numbers and the long-range coherence values. In figure [2] 
we compare the entropy calculated from the 6x6 density 
matrix estimate pcstice, (5) that is built by measuring sep- 
arately only the average well populations and the relative 
phases between the atoms in the 6 different sites 



Pest (a, /3) 



^(na}(n^)(exp[z((^a - h)]) 
(Nt) 



(5) 



where fia and c/jq, denote the atom number and phase op- 
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FIG. 3. Exact results for the total entropy St (solid lines) 
and its estimates (dashed lines) for different strengths of the 
barrier. The dashed black lines with crosses show the estimate 
derived from our approximate expression of the density matrix 
based on direct experimental observables ([5]); the red lines 
with circles show the estimate Sr from equation (|A20[I based 
on the Renyi entropies . Initial temperature T = 4nK for 
the ramping speeds (from left to right) 15.6 Hz/ms and 17.2 
Hz/ms. 



erators in the site a, respectively, and Nt represents the 
total atom number in the six sites. The good agreement 
between the estimated entropy and the one based on the 
full basis state representation can be explained by neg- 
ligible atom number pair correlations between different 
lattice sites {fiahf}) ~ {ha){np) and correlations between 
the phases and atom numbers {^Jfiafip exp[i(0Q— 0^)]) ~ 

(\/ nafip) {exp[i{(l)a — '/'/s)])- Similar correlations between 
the atoms in vibrational states within the same site only 
weakly affect the entropy. In addition, the loss of phase 
coherence between the vibrational levels k^l within the 
same site is small {eyi])[i{(j)ak—4'ai)\) ^ 0.99. The entropy 
approximation ^ is better at low temperatures owing 
to the weaker effect of thermal fluctuations on intra-site 
correlations between the atoms in different vibrational 
states. Analogously, stronger quantum fluctuations at 
stronger nonlinearities can lead to larger deviations from 
the exact result. Well-known estimates of von Neumann 
entropy are based on combinations of Renyi entropies [2^ 
(see,). It is therefore interesting to compare our estimate 
to the Renyi entropy estimates. We show in figure[3]how 
our entropy estimate based on experimental observables 
provides for this system a more accurate approximation 
of the entropy than the estimate Sr based on the Renyi 
entropies. 

Injecting laser barriers affects the achievable spin 
squeezing between the central sites after turning up of 
the lattice. We define the relative atom number squeez- 
ing between the atoms in sites a and /3 by 



= [A {ha - n/j)] 



{hg) + {hp) 
4:{ha){hp) 



(6) 



where A (ho.— hp) denotes the relative atom number 
fluctuations. The spin squeezing of the atoms. 
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not only depends on but also on the rela- 

tive phase coherence (cos(0^ ~ 4>a)) between the atoms. 
The proposed long-range phase coherence measurement 
scheme allows the detection of spin squeezing also be- 
tween non-nearest-neighbour sites and other spatially 
separated regions. In quantum-enhanced metrology a 
high-precision quantum interferometer can overcome the 
standard quantum limit of classical interferometers, pro- 
vided that Cs,(q,/3) < 1 [l^]- The same condition also 
implies quantum many-body entanglement in the sys- 
tem "3^. 

Results for the spin squeezing between the atoms in the 
two adjacent central wells (sites 3 and 4 of figure[T|) at the 
end of the ramping for two different ramping speeds are 
shown in figure Ul^a) and (b) at the experimentally rele- 
vant initial temperature T = 4nK [15]j. At intermediate 
scissor strengths Sf,, when the nonadiabatic injection does 
not significantly perturb the system [cf. figure [^Ja)], we 
find slightly stronger spin (as well as the relative atom 
number) squeezing than in the system where no laser 
potential was applied. The spin squeezing between the 
atoms in non-nearest-neighbour sites (sites 2 and 5 of 
figure [ll is shown in figure lUJc) and (d) . The system 
exhibits spin squeezing and quantum many-body entan- 
glement between spatially separated regions specified by 
the distant sites. Weak excitations of the system due to 
lattice and barrier ramping affect the squeezing as can be 
seen from the differences between the two ramping speed 
cases. At stronger scissor strengths the barriers perturb 
the system, resulting in notably stronger dynamics of the 
squeezing. We also find that in the region of improved 
spin squeezing the relative phase coherence is not no- 
tably affected. Our analysis provides a proof-of-principle 
demonstration that adjusting the coupling between the 
inner and outermost sites could potentially lead to tech- 
nologies for improving atomic spin squeezing. 

Spin squeezing flags the presence of entanglement in 
the system. As an alternative signature of entangle- 
ment one can use entropic inequalities in a bipartitite 



system 3ll433|. The purity Tr(p^(.^j) of the reduced den- 



sity matrices of the left (right) part of the system, com- 
pared to the purity of the combined system Tr[(pi 
signals entanglement if Tr(p^j^p < Tr[(p; ® Pr)^]- Even 
though we find spin squeezing in our system, the purity 
based entanglement inequality is not violated (seCp). Re- 
markably, however, the purity estimates derived from ([5]) 
agree well (figure[S]) with the exact results (within 1 % for 
Sh < 20 and within 6 % for all values of s^), thus poten- 
tially providing a direct experimental detection method 
for the presence of entanglement. 



^ Note that the actual experimental temperature in 3D trap may 
be higher than the one corresponding to experimental findings 
in ID simulations. 
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FIG. 4. Top row: Spin squeezing between the atoms in the 
two central lattice sites (sites 3 and 4 in figure [T| for difTerent 
strengths of the scissors potential, with the ramping speeds 
(from left to right) 15.6 Hz/ms and 17.2 Hz/ms. Bottom row: 
Spin squeezing between the atoms in non-nearest-neighbour 
sites (sites 2 and 5 in figure[l]), with the ramping speeds (from 
left to right) 15.6 Hz/ms and 17.2 Hz/ms. The dashed line 
shows the value obtained with no laser beam potential and 
the grey shaded area represents its uncertainty. The lattice 
is slowly ramped up to 72Er, after the barriers are rapidly 
introduced in thermal equilibrium at the lattice height 2AEr 
and T = 4nK. 



In conclusion, we have shown how the single-particle 
entropy can be measured and engineered in an optical 
lattice in the large filling limit. The measurement can 
be achieved by matter- wave homodyne detection of long- 
range phase coherence using a reference condensate. The 
proposed scheme provides a good approximation of the 
density matrix such that it might be used to detect entan- 
glement via entropic inequalities. We provided a proof- 
of-principle demonstration that controlling tunneling in 
selected locations of the lattice may find quantum tech- 
nological applications in manipulating metrologically im- 
portant spin squeezing. Improvement in spin squeezing 
has practical implications since in the experiments it is 
typically limited by thermal fluctuations 0. We may 
also envisage a procedure where the atomic spin squeez- 
ing is improved during every injection cycle of the laser 
barriers that is followed by an adiabatic adjustment of 
the trapping potential. Iterating the process could then 
potentially lead to progressively stronger spin squeezing. 



APPENDIX 
Stochastic phase-space method 
Dynamics 

In order to simulate the experimental conditions, quan- 
tum and thermal fluctuations are included via a classical 
stochastic fleld description based on TWA [IB-iil. TWA 
has provided a useful methodology for analyzing rela- 
tive atom number correlations between different lattice 
sites and, e.g., between different solitons in interfer- 



ometric applications 3J| . The dynamics is unraveled into 
stochastic trajectories where the initial state of the quan- 
tum operator iplxjt = 0) is represented by an ensemble 
of classical flelds ipw{x,t = 0), sampled according to its 
Wigner distribution that is determined by the temper- 
ature, the atom statistics and the multimode nature of 
the system. The initial state ^j{x,t = 0) is expanded in 
terms of the BEC ground state and the excited states as 

-(/'(x, t = Q) = V'oao + ^ |^Uj(a;)dj — w*(x)d] . (A8) 
i>o 

Here olj represent the quasiparticle annihilation opera- 
tors and the quasiparticle mode functions Uj (x) , Vj (x) 
are calculated by solving self-consistently the coupled 
Hartree-Fock-Bogoliubov equations 27|, |28| for the con- 
densate and non-condensate populations. The quasipar- 
ticle operators are than replaced by stochastic complex 
variables aj , a* , obtained by sampling the correspond- 
ing Wigner distribution. Each individual stochastic re- 
alization is dynamically evolved according to the Gross- 
Pitaevskii equation and represents a possible outcome of 
an individual experimental run. Ensemble averages cal- 
culated from the TWA numerical results give a statistical 
description of the dynamics of the system. 



Analysis of correlations 

The multi-mode nature of the system and phonon- 
phonon interactions were shown to be important for the 
evaluation of the atom number fluctuations 151. Multi- 



mode effects are included by projecting the stochastic 
field il^w onto a mode function basis formed by energy 
eigenfunctions of the lattice sites for a given number 
of energy bands. The projection technique allows to 
transform the symmetrically-ordered expectation values 
of stochastic representations of quantum operators in the 
Wigner distribution to normally-ordered expectation val- 
The projected quantities are used in 



ues 



calculation of the atom number statistics and the phase 
coherence of the system. This also provides a model 
for the single particle entropy, evaluated from the single- 
particle density matrix in the same mode function basis. 
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We define an eigenmode basis for each site given by 
the mode functions iprjj{x) with i] = \,...,Ny^,j = 
1, . . . ,Nm- Here the first index runs over all the sites, 
with = 6 in our case, and the second index runs 
over all the vibrational state mode functions in individ- 
ual sites. We denote the stochastic amplitude for the 
atoms in the jth vibrational mode of site 77 as a,; which 
can be numerically obtained from the projection of the 
stochastic field tpw{x^t) as 



(A9) 

The macroscopic phase for the atoms in each site may be 
calculated by averaging over the vibrational states 

<t>n{t) -arg / V a^j(i)93^.j (x) dx . (AlO) 

Jj)«'well J 

The relative phase coherence between the atoms in the 
sites r] and fj, can then be obtained from 



^'t'rip. = (cos {4>ri - 4'p.))w , 



(All) 



where the subscript W denotes the Wigner expectation 
value over many realizations. The projected amplitudes 
may be used to calculate the various normally ordered 
quantum expectation values. The site populations read 

(^'7) = = J2 [ialjar,,j)w - 1/2], (A12) 



where the summation is over all the vibrational modes in 
the site ry. The contribution —1/2 in the last term is a 
result of the symmetrical ordering of Wigner expectation 
values. The on-site atom number fluctuations for site rj 
are similarly given by 

I )w{\an,i\ )w{\arj,k\ )w 



-<5,;,fc/4] 



(A13) 



whereas the relative atom number fluctuations between 
two sites r] and fi are obtained from 

l2 



At 



[A {flri - flf,)] 



(A14) 



-(I 



i,k 



I /I 1^ I l^\ 

\0,ti,i\ )w{\0.7j,k\ — \0,ti,k\ )W — 



The single-particle density matrix p can be evaluated 
from the projected amplitudes and it is given by the en- 
semble average of the stochastic realizations pk where 



(A15) 



The matrix elements p^^^ gi are given by 



{Nt) 



(A16) 



Here (-/Vt) denotes the atom number in the studied six 
central sites. The diagonal elements are determined by 
the mode populations and the off-diagonal elements con- 
tain information on the phase coherence. 



Entropy 

We first calculate the single-particle (generalized von 
Neumann [37[) entropy S of equation ^ numerically in 
the projected single-particle basis and then propose an 
estimate of the entropy based only on direct experimen- 
tal observables 0, [l^- After injecting the barriers the 
entropy per atom Si can be evaluated in the new reduced 
lattice system of the four sentral sites. 



Si = 



(A17) 



where the prime in the sum indicates a summation over 
the wells that are inside the scissors and S{x) the entropy 
density. We compare this with the entropy per atom St 
of the full system 



St — 



J2,jIr,^^^cllSix)dx 



(A18) 



where the summation is now over all the sites. 

For comparison we also calculate an estimate of the 
von Neumann entropy provided by a combination of the 
Reiwi entropies. The Renyi entropy of order n is given 

by m 



1 



1-n 



log(Tr(/)")), n>2. 



(A19) 



In the limit n 1, it coincides with the von Neumann 
entropy. The Renyi entropy can be used to approximate 
the von Neumann entropy [39| . Specifically, an estimate 
Sr of the von Neumann entropy in terms of the Renyi 
entropy can be given as 

Sr = ^Su + ^max [^^23, ^12] , (A20) 

where Su, Sd23, Sf2 are appropriately defined functions 
of 52 and ^3 (see equation (32) in [29|). We find that 
our entropy estimate that is based on experimental ob- 
servables provides a more accurate approximation of the 
entropy than the one based on the Renyi entropies, as 
shown in figure [3] 
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FIG. 5. Exact results (solid lines) and estimates (dashed lines) 
for the purity of system A (Trp^, red lines with crosses) and 
of the combined bipartite system C [Trpc, black lines with 
circles) for different strengths of the barrier for the ramping 
speeds (from left to right) 15.6 Hz/ms and 17.2 Hz/ms. The 
dashed lines show the estimates derived from our approximate 
expression of the density matrix based on direct experimental 
observables ([5]). Initial temperature T = 4nK. The relative 
sampling errors for the exact results are smaller than those 
shown in figure (2] 



Entanglement 



Constructing the density matrix has useful applica- 
tions, e.e^ for identifying bipartite entanglement in the 



system 



31, 3 



A simple test can be done by comparing 
the purity of the full system described hy pc = PA® Pb 
to that of a subsystem (described by the reduced density 
matrix Pa(b))- Separability of the two subsystems re- 
quires Tr^p^j.^^^ > Tr (/O^). Since there are only a few 
experimentally accessible tests for entanglement in many 
body systems, it is interesting to compare the purities ob- 
tained from the proposed estimate of the density matrix 
to the exact results. As an example, we define the left 
half of the lattice system as subsystem A, the right half 
as subsystem B. We calculate the purity of the reduced 
density matrix of the left-half subsystem A and compare 
it to the purity of the combined bipartite system. As 
shown in figure [5] this test does not detect entanglement, 
but remarkably the results calculated from our density 
matrix estimate ([S]) based on experimental observables 
provides a good approximation of the exact result: The 
purities agree within 1 % for low barrier heights (below 
Sb = 20) and within 6 % for all values of Sb- 
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